Biological  Report  90(10) 
October  1990 


Tests  of  a  Habitat  Suitability  Model 
for  Black-capped  Chickadees 


Biological  Report 

This  publication  series  of  the  Fish  and  Wildlife  Service  comprises  reports  on  the  results  of  research,  developments  in 
technology,  and  ecological  surveys  and  inventories  of  effects  of  land-use  changes  on  fishery  and  wildlife  resources. 
They  may  include  proceedings  of  workshops,  technical  conferences,  or  symposia;  and  interpretive  bibliographies. 


Copies  of  this  publication  may  be  obtained  from  the  Publications  Unit,  U.S.  Fish  and  Wildlife 
Service,  1849  C  Street,  N.W.,  Mail  Stop  130-ARLSQ,  Washington,  DC  20240,  or  may  be  purchased 
from  the  National  Technical  Information  Service  (NTIS),  5285  Port  Royal  Road,  Springfield,  VA 
22161. 


ISSN  0895-1926 


Biological  Report  90(10) 
October  1990 


Tests  of  a  Habitat  Suitability  Model 
for  Black-capped  Chickadees 


By 


Richard  L.  Schroeder 


U.S.  Department  of  the  Interior 
Fish  and  Wildlife  Service 
Washington,  D.C.  20240 


Contents 


Page 

Abstract  .  1 

Study  Area  and  Methods  .  2 

Volume  Test  .  2 

HSI  Model  Test  .  3 

Statistical  Analyses  .  3 

Results  .  4 

Volume  Test  .  4 

HSI  Model  Test  . ^  ^  5 

Discussion  .  5 

Volume  Test  .  5 

HSI  Model  Test  .  7 

Acknowledgments  .  g 

References  .  g 


ii 


Tests  of  a  Habitat  Suitability  Model  For  Black-capped  Chickadees 


by 


Richard  L.  Schroeder 

U.S.  Fish  and  Wildlife  Service 
National  Ecology  Research  Center 
4512  McMurray  Avenue 
Fort  Collins,  Colorado  80525 


ABSTRACT. — The  black-capped  chickadee  {Parus  atricapillus)  Habitat  Suitability  Index  (HSI) 
model  provides  a  quantitative  rating  of  the  capability  of  a  habitat  to  support  breeding,  based  on 
measures  related  to  food  and  nest  site  availability.  The  model  assumption  that  tree  canopy  volume 
can  be  predicted  from  measures  of  tree  height  and  canopy  closure  was  tested  using  data  from 
foliage  volume  studies  conducted  in  the  riparian  cottonwood  habitat  along  the  South  Platte  River  in 
Colorado.  Least  absolute  deviations  (LAD)  regression  showed  that  canopy  cover  and  overstory 
tree  height  yielded  volume  predictions  significantly  lower  than  volume  estimated  by  more  direct 
methods.  Revisions  to  these  model  relations  resulted  in  improved  predictions  of  foliage  volume. 
The  relation  between  the  HSI  and  estimates  of  black-capped  chickadee  population  densities  was 
examined  using  LAD  regression  for  both  the  original  model  and  the  model  with  the  foliage  volume 
revisions.  Residuals  from  these  models  were  compared  to  residuals  from  both  a  zero  slope  model  and 
an  ideal  model.  The  fit  model  for  the  original  HSI  differed  significantly  from  the  ideal  model,  whereas 
the  fit  model  for  the  revised  HSI  did  not  differ  significantly  from  the  ideal  model.  However,  both  the  fit 
model  for  the  original  HSI  and  the  fit  model  for  the  revised  HSI  did  not  differ  significantly  from  a 
model  with  a  zero  slope.  Although  further  testing  of  the  revised  model  is  needed,  its  use  is  recom¬ 
mended  for  more  realistic  estimates  of  tree  canopy  volume  and  habitat  suitability. 


The  U.S,  Fish  and  Wildlife  Service  (Service)  has  pub¬ 
lished  Habitat  Suitability  Index  (HSI)  models  (Scham- 
berger  et  al.  1982)  for  more  than  150  species  of  fish  and 
wildlife.  These  models  were  developed  primarily  from 
literature  sources  and  expert  reviews  and  provide  a  quan¬ 
titative  method  to  predict  habitat  suitability  based  on  vege¬ 
tative  and  physical  characteristics  of  habitats.  Whereas 
earlier  Service  efforts  were  directed  at  model  develop¬ 
ment,  the  emphasis  in  recent  years  has  shifted  to  tests  of 
published  model  hypotheses  (Schamberger  and  O’Neil 
1986),  I  present  results  from  tests  of  hypotheses  in  the 
black-capped  chickadee  {Parus  atricapillus)  HSI  model 
(Schroeder  1983),  as  well  as  an  empirical  test  of  the  HSI 
model  itself. 

The  black-capped  chickadee  HSI  model  provides  a 
quantitative  rating  of  the  capability  of  a  habitat  to  support 
breeding,  based  on  measures  related  to  food  and  nest  site 
availability  (Fig.  1;  Table  1).  The  food  portion  of  the 
habitat  model  was  developed  from  data  presented  in  Stur- 
man  (1968)  and  assumes  that  black-capped  chickadee 
population  densities  are  related  to  arthropod  abundance 
and  that  this  is  positively  correlated  with  tree  foliage  vol¬ 
ume,  Sturman  (1968)  found  that  canopy  volume  of  trees 
was  a  strong  predictor  of  black-capped  chickadee  abun¬ 


dance  (r^  =  0.796).  However,  measures  of  tree  canopy 
volume  as  presented  in  Sturman  (1968)  are  very  time 
consuming  and  may  not  be  suitable  for  many  users  of  the 
HSI  model.  An  alternative  method  is  presented  in  the  HSI 
model,  based  on  the  assumption  that  measures  of  tree 
heiGht  and  tree  canopy  closure  can  be  used  in  place  of 
Sturman ’s  volume  measurements  to  assess  the  suitability 
of  the  food  resource  for  black-capped  chickadees. 

The  nest  site  availability  portion  of  the  HSI  model  as¬ 
sumes  that  reproduction,  and  hence  black-capped  chick¬ 
adee  numbers,  may  be  limited  by  the  density  of  snags  that 
are  10  to  25  cm  diameter  at  breast  height  (dbh).  This  part  of 
the  model  incorporates  estimates  of  the  number  of  suitable 
snags  per  unit  area  that  are  assumed  to  be  required  by 
breeding  black-capped  chickadees. 

The  black-capped  chickadee  HSI  model  was  also  tested 
in  New  York  by  Bayer  and  Porter  (1988).  They  compared 
population  densities  and  HSI  model  values  for  28  sites 
using  both  continuous  and  discrete  (>0.5,  >0-0.5,  and  0) 
HSI  values.  The  model  was  successful  in  predicting  popu¬ 
lation  abundance  at  the  discrete  level  {P  <  0.05)  but  not  at 
the  continuous  level. 

The  objectives  of  this  study  were  ( 1 )  to  test  the  hypothe¬ 
sis  that  tree  canopy  volume  can  be  predicted  from  meas- 
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Fig.  1.  Habitat  variable  and  suitability  index  graphs  for  the 
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black-capped  chickadee  HSI  model  (Schroeder  1983). 


ures  of  tree  height  and  canopy  closure,  and  (2)  to  test  the 
relation  between  the  HSI  and  estimates  of  black-capped 
chickadee  densities  during  the  breeding  season. 

Study  Area  and  Methods 

Foliage  volume  studies  were  conducted  in  plains  cot¬ 
tonwood  (Populus  sargentii)  bottomland  along  the  South 
Platte  River  near  Crook,  Logan  County,  Colorado.  Black- 
capped  chickadee  population  data  were  obtained  from 
previous  studies  in  the  same  area  conducted  by  Sedgwick 
and  Knopf  (1986;  J.  A.  Sedgwick  and  F.  L.  Knopf,  U.S. 
Fish  and  Wildlife  Service,  Fort  Collins,  Colorado,  un¬ 
published  data).  Elevations  along  the  South  Platte  River 
ranged  from  1 , 1 1 6  to  1 , 149  m,  and  plains  cottonwood  was 
the  dominant  overstory  tree  species.  For  further  descrip¬ 


tions  of  the  riparian  vegetative  community  refer  to  Sedg¬ 
wick  and  Knopf  (1986). 

Volume  Test 

A  major  assumption  in  the  HSI  model  is  that  the  food 
suitability  index  (SI)  computed  through  measuring  tree 
height  and  canopy  closure  is  equivalent  to  the  food  SI 


Table  1 .  Life  requisite  and  HSI  determination  methods 
for  the  black-capped  chickadee  (Parus  atricapillusj  HSI 
model  (Schroeder  1983). 


Model  output 

Determination  method 

Food  suitability  index  (SI) 

(SIVl  XSIV2)'/2  or  (SIV3) 

Nest  SI 

SIV4 

HSI 

Lower  of  food  SI  or  nest  SI 

Black-capped  Chickadees 
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value  computed  for  tree  canopy  volume  as  measured  using 
Sturman’s  (1968)  methods.  If  this  assumption  of  equiva¬ 
lence  is  true,  canopy  volume  can  be  predicted  from  tree 
height  and  canopy  closure  by  calculating  the  food  SI  for 
tree  height  and  canopy  closure  and  then  determining  the 
canopy  volume  associated  with  the  food  SI  from  the  graph 
of  SIV3  (tree  canopy  volume).  A  food  SI  of  1.0  produced 
through  measuring  tree  height  and  canopy  closure  should 
indicate  a  tree  canopy  volume  >10.2  m^/m^  of  ground 
surface.  The  null  hypothesis  for  the  tree  canopy  volume 
test  is  that  for  tree  canopy  volumes  <10.2  m^m^  ground 
surface  there  is  a  linear  ratio  of  1 0.2: 1  between  tree  canopy 
volume  predicted  using  the  methods  of  Sturman  (1968) 
and  the  food  SI  determined  through  the  use  of  tree  height 
and  tree  canopy  closure,  with  an  intercept  of  0. 

To  test  this  hypothesis,  detailed  measurements  were 
taken  of  each  tree  (minimum  height  of  5  m)  on  plots 
selected  to  provide  wide  variation  in  the  amount  of  foliage 
volume  estimated  to  be  present.  Eighteen  plots  were  se¬ 
lected  with  plot  sizes  of  0.04  ha  except  for  3  0. 1  -ha  plots  in 
scattered  mature  cottonwood  stands  and  2  0.01 -ha  plots  in 
dense  cottonwood  sapling  stands. 

Measurements  to  determine  the  volume  of  living  fo¬ 
liage  of  each  tree  were  taken  in  May  1988,  following  the 
methods  described  in  Sturman  ( 1 968).  Tree  height  meas¬ 
urements  were  taken  with  a  clinometer  at  a  measured 
distance  from  each  tree.  The  width  of  individual  tree  cano¬ 
pies  was  measured  by  vertically  projecting  the  inner  and 
outer  portion  of  the  living  foliage  and  marking  these  points 
on  the  ground  with  wire  flags. 

The  percent  of  tree  canopy  cover  within  each  plot  was 
measured  with  a  spherical  densitometer  (Lemmon  1956) 
at  20  points,  equally  spaced  along  2  diagonal  transects. 
The  densitometer  was  modified  to  reduce  the  overhead 
angle  viewed  in  the  mirror  by  using  only  the  four  center 
squares.  This  technique  ensured  that  canopy  cover  read¬ 
ings  were  nearly  vertical  and  did  not  take  in  vegetation 
outside  the  plot.  A  hit  or  miss  was  recorded  on  each  of  four 
points  within  each  mirror  square. 

Measurements  of  tree  dbh  were  taken  to  determine 
basal  area  for  each  plot.  Basal  area  measurements  were  not 
required  in  the  original  HSI  model.  However,  basal  area  is 
widely  used  and  was  measured  to  determine  the  relation 
between  it  and  tree  canopy  volume. 

HSI  Model  Test 

The  null  hypothesis  for  the  test  of  the  HSI  model  is  that 
regression  of  the  HSI  on  estimated  black-capped  chick¬ 
adee  population  densities  will  produce  a  model  with  a 
slope  significantly  different  from  zero  but  not  signifi¬ 
cantly  different  from  the  slope  of  a  model  predicting  max¬ 
imum  expected  population  densities.  Black-capped  chick¬ 
adee  breeding  population  data  from  previous  studies 
(Sedgwick  and  Knopf  1986  and  unpublished  data)  were 


used  to  test  this  hypothesis.  These  data  were  collected  on 
10  16-ha  plots  in  riparian  cottonwood  stands  in  1982  and 
1984-86,  using  the  variable  circular-plot  method  (Rey¬ 
nolds  et  al.  1980).  Black-capped  chickadee  population 
density  estimates  were  determined  by  using  program 
TRANSECT  (Burnham  et  al.  1980).  Tree  heights  were 
measured  with  a  clinometer.  Data  on  snag  density,  as 
defined  in  the  HSI  model,  were  not  collected;  however, 
the  number  of  trees  with  cavities  (any  opening  3  to  12  cm 
in  diameter)  was  recorded.  Based  on  the  scarcity  of  snags 
on  the  study  site  (0.66  snags  per  hectare),  and  the  fact 
that  most  cavities  (95.8%)  were  in  live  trees,  it  seems  that 
the  number  of  trees  with  cavities  provides  a  reasonable 
substitute  for  the  snag  density  data  required  in  the  HSI 
model.  Additional  descriptions  of  methods  for  population 
estimation  and  vegetation  measurement  are  presented  in 
Sedgwick  and  Knopf  (1986).  I  obtained  measurements 
of  tree  canopy  closure  from  low-level  aerial  photographs 
of  each  of  the  10  plots.  Individual  canopies  of  cottonwood 
trees  were  outlined  on  an  acetate  overlay  of  the  photo¬ 
graphs,  and  digitized.  Tree  canopy  cover  was  computed  by 
comparing  the  total  area  of  the  tree  canopies  to  the  area 
of  each  plot. 

Statistical  Analyses 

I  used  least  absolute  deviations  (LAD)  regression  (Bar- 
rodale  and  Roberts  1974)  to  examine  linear  relation  be¬ 
tween  (1)  tree  canopy  volume  predicted  from  tree  height 
and  canopy  closure,  and  canopy  volume  estimated  follow¬ 
ing  Sturman  (1968),  and  (2)  black-capped  chickadee  den¬ 
sities  and  HSI  values.  Least  absolute  deviation  regression 
is  more  resistant  to  the  influence  of  outlying  data  values 
than  least  squares  (LS)  regression  (Wilson  1978;  Narula 
and  Wellington  1982),  and  LAD  regression  estimators 
have  smaller  variances  than  LS  estimators  for  data  that 
deviate  from  normality  and  homogeneous  variance,  espe¬ 
cially  with  small  sample  sizes  (Wilson  1978;  Dielman  and 
Pfaffenberger  1982).  Plots  of  these  data  indicated  nonnor¬ 
mality  and  that  heterogeneous  variances  were  present. 
Least  squares  regression  lines  and  values  are  given  for 
comparative  purposes. 

Statistical  inferences  about  proposed  and  fit  regression 
models  (Ho:  pfit  =  pproposed;  Ha:  Pfit  ^  Pproposed)  were 
made  by  using  a  distribution-free  permutation  test  for 
matched  pairs  (PTMP;  Mielke  and  Berry  1982)  to  com¬ 
pare  absolute  values  of  the  LAD  residuals  for  alternative 
models.  Euclidean  distances  were  used  in  the  PTMP 
analyses  because  they  are  consistent  with  the  geometry  of 
the  LAD  regressions  (Mielke  1986),  The  combination  of 
LAD  regression  and  distribution-free  tests  of  the  LAD 
residuals  is  a  more  desirable  approach  than  transforming 
variables  to  meet  LS  assumptions  because  it  avoids  arbi¬ 
trarily  rescaling  data  into  a  geometry  that  is  inconsistent 
with  the  original  measurement  scale. 
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Results 

Volume  Test 

Least  absolute  deviation  regression  lines  forced 
through  an  intercept  of  0  (because  there  can  be  no  volume 
if  there  are  no  trees)  produced  the  following  fit  model 
(Fig.  2): 

Y^3.21(X), 

where  Y  -  volume  (following  Sturman’s  methods), 

X  =  food  SI,  computed  as  (SIVl  X  SIV2)‘/2. 

According  to  the  HSI  model,  using  Sturman’s  (1968)  data, 
the  proposed  model  should  be  of  the  form  Y  =  10.2(X). 
Residuals  differed  between  the  fit  and  proposed  models 
(P  -  0.003),  indicating  significant  differences  between 
slopes  of  the  fit  and  the  proposed  models.  Plots  of  these 
data  (Fig.  2),  along  with  least  squares  regression  analy¬ 
sis  (Y  =  3.9 1[X],  /*2  =  0.70),  indicate  that  a  positive  re¬ 
lation  exists  between  the  model  output  and  measured 
foliage  volume. 

The  plots  show  that  the  HSI  model  assigns  high  index 
values  to  plots  that  have  relatively  low  levels  of  foliage 
volume.  The  following  changes  were  made  in  the  model  to 
improve  the  relation  (Fig.  3).  First,  the  SI  graph  for  tree 
canopy  closure  was  modified  to  obtain  maximum  Si’s  at 
>70%  canopy  closure.  The  original  model  (Fig.  1)  pro¬ 
posed  a  maximum  SI  for  this  variable  between  50  and  75% 
canopy  closure,  with  decreasing  values  above  75%.  These 
data  indicate  that  this  assumption  was  not  correct  for  the 
habitat  type  studied.  The  modified  relation  seems  to  more 
accurately  depict  the  effects  of  changes  in  canopy  closure 
on  foliage  volume. 


COTTONWOOD 


PROPOSED  MODEL 
Y=10.2(X) 


LAD  REGRESSION 
Y  =  3,21  (X) 


FOOD  Si  (square  root  of  SIV1  x  SIV2) 

Fig.  2.  Plots  of  original  food  suitability  indices  for  cottonwood 
habitats  versus  tree  canopy  volume  measured  following  Stur- 
man  (1968). 


Second,  the  original  variable  “average  height  of  over¬ 
story  trees”  considered  only  the  height  of  overstory  trees 
(defined  as  trees  >80%  of  the  height  of  the  tallest  tree, 
where  trees  are  considered  to  be  any  woody  vegetation  >5 
m  in  height).  Thus,  the  indirect  estimates  of  canopy  vol¬ 
ume  in  stands  with  a  few  very  tall  trees  and  many  trees 
<80%  of  the  height  of  the  tallest  tree  would  be  overrated 
with  the  original  model  because  the  limited  volumes  of  the 
many  smaller  trees  would  not  be  averaged  with  the  taller 
trees.  Therefore,  the  model  was  revised  to  consider  the 
average  height  of  all  trees.  The  original  model  attained 
maximum  SI  values  at  tree  heights  >15  m.  These  data 
indicate  that  volume  continues  to  increase  as  tree  heights 
increase  above  15  m.  Therefore,  the  SI  relation  was  mod¬ 
ified  to  provide  increasing  values  as  tree  heights  increased 
up  to  25  m.  In  addition,  the  original  model  provided  SI 
values  for  woody  vegetation  of  any  height  >0  m.  Thus,  a 
stand  with  an  average  tree  height  of  5  m  would  receive  an 


SIVl 


SIV2 


AVERAGE  HEIGHT  OF  TREES  (m) 


Fig.  3.  Revised  suitability  index  graphs  for  SIVl  and  SIV2. 
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SI  of  0.33.  These  stands  would  not  provide  suitable  breed¬ 
ing  habitat  for  the  black-capped  chickadee,  which  is  a 
forest-dwelling  bird.  Therefore,  the  SI  relation  was  mod¬ 
ified  to  provide  SI  values  of  0  for  stands  with  average  tree 
heights  of  <5  m. 

The  last  change  made  in  the  model  was  to  eliminate  the 
square  root  function  from  the  food  SI  calculation.  The 
volume  of  a  three-dimensional  object  is  equal  to  the  prod¬ 
uct  of  the  area  of  the  base  and  the  height.  For  a  forest 
stand,  a  measure  of  tree  canopy  closure  is  comparable  to 
the  area  occupied  by  tree  canopies,  and  a  measure  of 
average  tree  height  considers  the  vertical  height.  Thus,  it  is 
more  logical  to  compute  the  food  SI  by  using  the  product 
of  SIVl  (canopy  area)  and  SIV2  (height),  rather  than  the 
geometric  mean. 

These  model  revisions  had  a  marked  effect  on  the  abil¬ 
ity  of  the  model  to  predict  tree  foliage  volume.  Least 
absolute  deviation  regression  analysis  of  the  revised  food 
SI  relations,  forcing  a  0  intercept,  produced  the  following 
fit  model  (Fig.  4): 

Y=  10.1 1(X), 

where  Y  =  volume  (following  Sturman’s  methods), 

X=  food  SI,  computed  as  (Rev.  SIVl  X  Rev. 

SIV2). 

Residuals  from  the  fit  and  proposed  (Y  =  10.2[X])  models 
did  not  differ  (P  =  0.61),  indicating  no  significant  differ¬ 
ences  between  slopes  of  the  fit  and  proposed  models.  The 
least  squares  regression  fit,  again  with  intercept  forced  to 
0,  was  Y  =  10.39(X)  =  0.877). 

Data  on  basal  area  for  each  plot  were  also  regressed 
against  volume  as  computed  using  Sturman’s  (1968) 
methods,  and  it  seems  that  basal  area  is  a  good  predictor  of 
foliage  volume.  Least  absolute  deviation  regression  analy¬ 


sis,  with  the  intercept  forced  to  0,  produced  the  following 
fit  model: 

Y  =  0.138(X), 

where  Y  =  volume  (following  Sturman’s  methods), 

X  =  basal  area  from  plot  data  (m^/ha). 

Residuals  from  the  fit  models  and  models  with  a  slope  of 
zero  were  different  (P  =  0.05),  indicating  the  slope  of  the 
fit  model  was  significantly  different  than  zero. 

HSI  Model  Test 

The  test  of  the  volume  relation  noted  previously  only 
considers  the  model’s  ability  to  predict  tree  canopy  vol¬ 
ume  based  on  measures  of  tree  height  and  canopy  closure 
or  basal  area.  A  significant  additional  question  is  the  rela¬ 
tion  of  tree  canopy  volume  to  measures  of  black-capped 
chickadee  abundance.  The  hypothesis  in  the  HSI  model  is 
that  habitats  will  provide  increasing  potential  black- 
capped  chickadee  population  densities  as  the  HSI  in¬ 
creases  from  0  to  1 . 

Both  the  original  and  the  revised  food  SI  relations  were 
used  to  test  the  relation  between  the  HSI  and  black-capped 
chickadee  population  densities  from  the  10  16-ha  riparian 
cottonwood  plots.  The  model  states  that  the  HSI  equals  the 
lower  of  either  the  food  or  nest  SI  values.  Population 
densities,  original  and  revised  food  Si’s,  nest  Si’s,  and 
original  and  revised  HSI  values  for  the  10  plots  are  pro¬ 
vided  in  Table  2.  Least  absolute  deviation  regressions  of 
population  densities  on  the  original  and  revised  HSI  val¬ 
ues  produced  the  following  fit  models: 

Original  HSI  Model: 


FOOD  SI  (SIVl  X  SIV2) 

Fig.  4.  Plot  of  revised  food  suitability  indices  for  cottonwood 
habitats  versus  tree  canopy  volume  measured  following  Stur- 
man  (1968). 


Y  =  0.098 +  0.416(X), 

where  Y  =  population  density  of  black-capped  chick¬ 
adees  (both  sexes)  per  hectare, 

X  =  original  HSI. 

Revised  HSI  Model: 

Y  =  0.096  +  2.344(X), 

where  Y  =  population  density  of  black-capped  chick¬ 
adees  (both  sexes)  per  hectare, 

X  =  revised  HSI. 

Residuals  of  the  fit  models  and  a  model  with  a  zero  slope 
did  not  differ  (original  HSI  model,  P  =  0.20;  revised  HSI 
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Table  2.  Black-capped  chickadee  fParus  atricapillusj  densities,  food  and  nest  SI  values,  and  MSI  values  for  the  original  and  revised 

models. 


Plot 

number 

Chickadees 
per  hectare 

Original 
food  SI 

Revised 
food  SI 

Nest 

SI 

Original 

HSI 

Revised 

HSI 

1 

0.23 

0.62 

0.097 

0.88 

0.62 

0.097 

2 

0.39 

0.70 

0.084 

0.81 

0.70 

0.084 

3 

0.46 

0.62 

0.080 

0.63 

0.62 

0.080 

4 

0.37 

0.81 

0.102 

0.57 

0.57 

0.102 

5 

0.30 

0.58 

0.087 

0.39 

0.39 

0.087 

6 

0.30 

0.68 

0.113 

0.41 

0.41 

0.113 

7 

0.29 

0.60 

0.071 

0.58 

0.58 

0.071 

8 

0.10 

0.58 

0.048 

0.19 

0.19 

0.048 

9 

0.14 

0.47 

0.052 

0.10 

0.10 

0.052 

10 

0.15 

0.49 

0.021 

0.31 

0.31 

0.021 

model,  P  =  0.34),  indicating  that  neither  of  the  fit  models 
differed  significantly  from  a  zero  slope. 

An  important  aspect  of  testing  HSI  models  is  to  com¬ 
pare  the  model  output  to  an  estimate  of  maximum  ex¬ 
pected  performance  of  the  species.  The  black-capped 
chickadee  HSI  model  was  developed  to  encompass  the 
entire  range  of  the  species.  To  determine  the  highest  popu¬ 
lation  densities  that  have  been  encountered  in  previous 
studies,  data  from  more  than  350  Breeding  Bird  Census 
(BBC)  reports  (T.  Engstrom,  Cornell  Laboratory  of  Or¬ 
nithology,  Ithaca,  New  York,  personal  communication) 
were  analyzed.  These  data  indicated  that  the  average  den¬ 
sity  of  the  highest  five  BBC ’s  for  black-capped  chickadees 
was  1 .4  males  per  hectare  (excluding  one  small  orchard 
with  1  breeding  male  and  a  corresponding  density  of  2.5 
males  per  hectare).  Assuming  that  actual  density  of  both 
sexes  would  be  double  this  figure  (2.8  chickadees  per 
hectare),  a  proposed  model  for  an  ideal  HSI  would  have 
the  form 

Y  =  2.8(X), 

where  Y  =  population  density  of  black-capped  chick¬ 
adees  (both  sexes)  per  hectare, 

X  =  HSI. 

Residuals  from  the  fit  model  for  the  original  HSI  and  the 
proposed  ideal  model  were  different  (P  =  0.002),  indicat¬ 
ing  a  significant  difference  between  the  slopes  of  these 
models.  This  difference  can  be  seen  because  the  fit  model 
for  the  original  HSI  yields  a  predicted  population  density 
of  0.52/ha  when  the  HSI  is  at  its  maximum  value  of  I .  Re¬ 
siduals  from  the  fit  model  for  the  revised  HSI  and  the  pro¬ 
posed  ideal  model  did  not  differ  (P  =  0.50),  indicating  no 
significant  difference  between  the  slopes  of  these  models. 
The  fit  model  for  the  revised  HSI  predicts  a  population 
density  of  2.44/ha  when  the  HSI  is  at  1.  Thus,  the  revised 
HSI  model  more  closely  approximates  maximum  ex¬ 
pected  population  densities  (2.8/ha)  when  the  HSI  is  at 


its  maximum.  The  LS  regression  fit  for  the  revised  HSI 
model  was  Y  =  0.063  +  2.772(X)  (r^  =  0.425),  yielding  a 
maximum  estimated  population  density  of  2.84yTia  when 
the  HSI  is  1. 


Discussion 

Volume  Test 

Proper  evaluation  of  HSI  models  involves  assessing  the 
performance  of  both  the  HSI  and  other  testable  assump¬ 
tions  within  the  model.  Schamberger  and  O’Neil  (1986) 
noted  that  most  studies  tested  the  HSI  against  some  meas¬ 
ure  of  population  status  (Clark  and  Lewis  1983;  Cook  and 
Irwin  1985;  Lancia  and  Adams  1985).  Here,  I  have  as¬ 
sessed  the  performance  of  both  a  model  assumption  and 
the  HSI  value. 

The  foliage  volume  portion  of  the  tests  indicates  that  the 
original  assumed  relation  of  canopy  cover  and  overstory 
tree  height  to  foliage  volume,  although  showing  a  strong 
linear  relation,  yielded  predictions  significantly  lower 
than  the  volume  estimated  by  Sturman’s  (1968)  methods. 
Revisions  to  these  model  relations  resulted  in  improved 
predictions  of  foliage  volume. 

Habitat  Suitability  Index  model  testing  can  lead  to  new 
information  if  the  test  is  not  limited  to  hypotheses  con¬ 
cerning  only  the  variables  found  in  the  HSI  model.  The 
measurement  of  tree  dbh  was  added  to  this  study  to  assess 
the  possibility  that  basal  area  could  be  used  to  predict  tree 
foliage  volume.  Basal  area  data  are  widely  available  and 
are  often  used  in  models  to  predict  future  forest  condi¬ 
tions.  Thus,  it  would  be  a  useful  variable  for  biologists 
conducting  forest  habitat  studies.  Data  for  the  cottonwood 
habitats  indicate  that  basal  area  is  an  excellent  predictor  of 
foliage  volume.  Based  on  these  results,  users  of  the  black- 
capped  chickadee  HSI  model  have  three  reasonable  op¬ 
tions  to  estimate  foliage  volume:  Sturman’s  (1968) 
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methods,  the  revised  SI  graphs  and  combination  formula 
for  tree  height  and  canopy  closure,  and  the  predictive 
model  using  basal  area. 

Sturman  (1968)  provided  only  two  potential  tree  shapes 
to  be  used  in  calculating  foliage  volume  for  an  individual 
tree,  one  each  for  conifers  and  deciduous  trees.  Many 
cottonwood  trees  do  not  clearly  match  these  shapes.  Given 
the  irregular  shapes  of  many  trees,  I  made  visual  estimates 
and  adjustments  of  the  outer  and  inner  canopy  layers  to 
approximate  the  desired  shape.  Mawson  et  al.  (1976)  de¬ 
scribed  a  method  to  determine  tree  crown  volume  using  15 
tree  shapes.  Improved  estimates  of  volume  may  be  possi¬ 
ble  by  expanding  the  number  of  tree  shapes  to  match  more 
realistically  what  is  encountered  in  forested  habitats, 

HSI  Model  Test 

Results  of  the  test  of  the  HSI  as  a  predictor  of  black- 
capped  chickadee  population  density  are  not  as  easily 
interpreted  as  the  volume  test  results.  The  fit  model  for  the 
revised  HSI  did  not  differ  from  the  ideal  model,  whereas 
the  fit  model  for  the  original  HSI  differed  significantly 
from  the  ideal  model.  However,  both  the  fit  model  for  the 
original  HSI  and  the  fit  model  for  the  revised  HSI  did  not 
differ  significantly  from  a  model  with  a  zero  slope.  Given 
the  changes  made  in  the  revised  HSI  model  and  the  spe¬ 
cific  data  set  used  to  test  the  model,  several  possible  expla¬ 
nations  can  be  offered  for  these  results:  (1)  Whereas  the 
revised  model  seems  to  be  a  better  predictor  of  tree  canopy 
volume,  volume  may  be  an  inadequate  predictor  of  popu¬ 
lation  densities  in  cottonwood  sites.  Additional  factors 
may  need  to  be  added  to  the  model  to  improve  its  predic¬ 
tive  ability;  (2)  the  number  of  sites  (n  =  10)  was  inadequate 
to  effectively  test  the  HSI  model.  If  the  same  spread  as  in 
these  data  occurred  over  30  sites,  the  revised  fit  model 
would  differ  significantly  from  a  zero  slope  {P  =  0.047), 
and  thus  fail  to  reject  both  portions  of  the  null  hypothesis. 
In  this  same  situation  (if  n  =  30)  the  original  fit  model 
would  continue  to  be  significantly  different  from  the  pro¬ 
posed  ideal  model  {P  <  0.0001 ),  thus  rejecting  this  portion 
of  the  null  hypothesis;  and  (3)  the  10  sites  used  for  the 
black-capped  chickadee  test  had  too  narrow  a  range  of  tree 
canopy  cover  values  (1 1  to  33%)  and  average  height  of 
trees  (7.5  to  12.3  m)  to  provide  an  adequate  test  of  the 
hypothesis.  Without  further  study,  it  is  not  possible  to 
determine  which  of  these  possibilities  best  explains  the 
results  obtained  in  this  test. 

Only  a  single  measure  of  black-capped  chickadee  per¬ 
formance  (number  per  hectare)  was  used  in  this  test.  How¬ 
ever,  population  density  may  be  a  misleading  indicator  of 
habitat  quality,  as  noted  by  Van  Home  (1983).  Additional 
measures  of  performance,  such  as  nest  success  or  number 
of  young  fledged,  might  improve  the  test.  In  addition,  the 
HSI,  as  a  measure  of  key  habitat  variables,  considers  only 


a  subset  of  all  of  the  factors  that  may  influence  individuals 
and  resultant  populations.  Thus,  habitat  measurements 
may  indicate  ideal  conditions,  but  nonhabitat-related  fac¬ 
tors  (such  as  disease  or  severe  weather)  may  cause  popula¬ 
tion  declines.  The  HSI  is  best  viewed  as  an  index  of 
habitat-imposed  limitations  on  a  population.  For  example, 
an  HSI  of  1  would  indicate  no  habitat-related  limitations 
exist,  whereas  decreasing  HSI  values  would  indicate  in¬ 
creasingly  significant  habitat  limitations.  Actual  popula¬ 
tion  densities  should  not  exceed  the  levels  indicated  by  the 
HSI  values  but  may  range  below  the  HSI  predictions  and 
not  falsify  this  hypothesis.  Specific  tests  related  to  this 
aspect  of  the  HSI  have  been  suggested  (U.S.  Fish  and 
Wildlife  Service  1987,  1989),  but  the  sample  size  in  this 
study  was  too  small  to  conduct  the  tests. 

One  suggested  change  in  the  nest  portion  of  the  black- 
capped  chickadee  HSI  model  resulted  from  this  study.  The 
original  model  required  a  measure  of  the  number  of  snags 
from  10  to  25  cm  dbh  as  an  indication  of  nest  site  avail¬ 
ability.  Sedgwick  and  Knopf  (1986),  however,  noted  that 
snag  densities  alone  do  not  accurately  reflect  the  availabil¬ 
ity  of  suitable  nesting  substrates  in  riparian  cottonwood 
stands.  The  majority  of  potential  nest  sites  in  these  stands 
occurred  in  cavities  in  live  trees.  In  eastern  hardwood- 
coniferous  forests,  use  of  live  trees  is  much  less  common 
(of  44  nests,  41  were  in  snags  or  stubs  and  only  3  were  in 
live  trees;  Runde  and  Capen  1987).  The  best  overall  meas¬ 
ure  of  black-capped  chickadee  nest  site  availability  may 
be  the  combined  density  of  the  number  of  trees  with  >1 
cavity  (minimum  diameter  of  10  cm)  and  the  number  of 
snags  (minimum  dbh  of  10  cm).  Further  studies  are  re¬ 
quired  to  test  whether  regional  differences  in  use  of  live 
trees,  as  opposed  to  snags,  should  be  incorporated  into  the 
HSI  model. 

I  have  not  developed  alternatives  to  the  revised  black- 
capped  chickadee  HSI  model,  nor  have  I  attempted  to 
search  out  the  best  model  (i.e.,  best  statistical  fit)  from  the 
data.  Fitting  a  new  model  to  these  1 0  data  points  may  result 
in  a  product  with  no  better  predictive  ability  than  the 
current  revision.  I  recommend  that  additional  tests  be  per¬ 
formed  on  the  revised  HSI  model  before  further  changes 
are  incorporated.  In  the  interim,  I  recommend  that  the 
revised  HSI  model  be  used  by  field  biologists  and  wildlife 
managers.  The  revised  model  uses  the  SI  relations  for 
canopy  closure  and  tree  height  illustrated  in  Fig,  3,  the 
food  SI  computed  as  the  product  of  SI VI  and  SIV2,  and 
the  variable  for  nest  suitability  (SIV4)  measured  as  the 
combined  density  of  the  number  of  trees  with  >1  cavity 
(minimum  diameter  of  10  cm)  and  the  number  of  snags 
(minimum  dbh  of  10  cm).  Although  additional  testing  is 
necessary,  it  seems  that  the  revised  model  will  provide 
more  realistic  estimates  of  both  tree  canopy  volume  and 
black-capped  chickadee  habitat  suitability  than  the  origi¬ 
nal  version. 
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